Cosmic String Induced CMB Maps 
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We compute maps of CMB temperature fluctuations seeded by cosmic strings using high resolu- 
tion simulations of cosmic strings in a Friedmann-Robertson-Walker universe. We create full-sky, 
18° and 3° CMB maps, including the relevant string contribution at each resolution from before 
recombination to today. We extract the angular power spectrum from these maps, demonstrating 
the importance of recombination effects. We briefly discuss the probability density function of the 
pixel temperatures, their skewness, and kurtosis. 



I. INTRODUCTION 

Despite improving observational limits, interest in cos- 
mic strings has remained durable (for a review, see 
Strings are a generic phenomena in fundamental theories 
and they can emerge in macroscopic form in braneworld 
cosmologies, for example, at the end of inflation 0, 
They are also common to cosmologically viable super- 
symmetry grand unified theory models _4] • Stringent con- 
straints on strings are important, therefore, in restricting 
the latitude available for cosmological model building. 
The detection of cosmic strings would be a watershed for 
high energy theory. 

Despite the potential significance, the investigation of 
cosmic strings and their observational consequences faces 
many numerical and analytic challenges, not least in cre- 
ating accurate realizations of string imprints in the cos- 
mic microwave sky. In this paper, we take this study 
a step further forward by presenting full-sky and small- 
angle CMB maps of temperature fluctuations seeded by 
cosmic string networks using high resolution simulations 
in an Friedmann-Robertson-Walker expanding universe 
(with the longest dynamic range to date). This work in- 
cludes all the relevant recombination physics and can be 
used not only to determine the angular power spectrum of 
string CMB anisotropics but also the higher order corre- 
lators such as the bispectrum, trispectrum, and beyond. 

Current constraints on cosmic strings result from line- 
of-sight CMB power spectrum calculations sourced either 
by unequal-time correlators obtained from field theory 
string simulations or semianalytic models of Nambu 
strings [1, |^. Qualitatively these two approaches pro- 
duce consistent spectra, that is, without the strong co- 
herent acoustic peaks associated with inflation. However, 
quantitatively there is a mismatch between the two ap- 
proaches in both the shape of the primary peak and its 
amplitude, which differs by a factor of 2-3. This dispar- 
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ity arises primarily from a difference in string network 
densities, which has been discussed at some length else- 
where [l^ (see also [HI). Nevertheless, there is general 
agreement that the relative amplitude of string induced 
CMB fluctuations cannot exceed more than 10% of those 
arising from adiabatic inflationary perturbations 13, 18j . 
There have also been a number of studies going beyond 
the power spectrum through map making with cosmic 
strings [l3 - [l5{ in order to study the degree of Gaus- 
sianity of the resulting CMB signatures. However, this 
work has generally only included late-time gravitational 
effects, ignoring the recombination physics which makes 
an important contribution to the signal over a wide range 
of muhipoles I « 200-2000. 

The motivation for the present work, then, is twofold: 
first, to include all recombination effects in the string 
CMB maps, so that we can ultimately characterize their 
primary statistical properties, and secondly, to match 
the accuracy of future experiments such as Planck H, 
AMI [13] and QUIET [l!] which will impose consider- 
ably more stringent constraints on cosmic strings through 
improving precision, resolution, and added polarization 
information. 



II. COSMIC STRING NETWORK 
SIMULATIONS 

Cosmic string simulations were performed with the 
AUen-Shellard string network code [l^. We have used 
fixed comoving resolution together with an initial string 
resolution of 24 points per correlation length. Simula- 
tions which started in the matter era had a dynamic 
range of 7.5 in conformal time, i.e. rjf = 7.5rii, but in 
order to use only the simulation when the network has 
relaxed into a scaling regime, we ignored the first 4% of 
time steps, resulting in an effective dynamic range of 6 in 
conformal time. Simulations that started in the radiation 
era had an effective dynamic range of 5, after eliminat- 
ing the first 4% of time steps from a simulation with a 
dynamic range of 6. The energy- momentum tensor of 
the network was projected onto a grid of 256'^ points as 
described in The background cosmology used was 



2 



ACDM with the WMAP 5-year data best fit parame- 
ters '2H: ncDM =0.214, =0.044, = 0.742, and h 
= 0.719. 

We have used three string simulations which span the 
range from before equality to today. The epochs of 
each simulation are (1) from rj = rjo/T.b to rjo, (2) from 
V = '7o/45 to 770/6, and (3) from 77 — 770/2I6 to 770/86. 
Thus, these simulations together span a range of 180 in 
conformal time (note that since we ignore the beginning 
of each simulation, they do not overlap). 



III. EINSTEIN-BOLTZMANN EVOLUTION 

We use the Landriau-Shellard code [2^ to com- 
pute cosmological perturbations in Fourier space. Two 
changes have been made to this code since the algorithm 
was presented. 

The first modification concerns the scalar metric equa- 
tions employed: Instead of solving for h and , we now 
solve for h and = h — , which obeys the following 
equation: 

h- + 2-h- = -167rG(pE + Sp) + SnGiQ + 29^) (1) 
a 

The other modification concerns the inverse compu- 
tation of the evolution equations' fundamental matrices: 
Instead of LU factorization, we use singular value de- 
composition, for which we employ the freely available 
LAPACK routines. This has proved a more numerically 
stable method and enables a better treatment of near 
singular matrices (22j . especially around recombination. 



In practice, because fe""^ and e~'^ are zero before the 
start of recombination, we only output grids for all per- 
turbations from rj = 27yrec/3 to the end of simulation 3, 
which is 77 ~ 3rirec/'2- For simulations 1 and 2, we only 
output the grids for hij, because fe""^ is also zero after 
the end of recombination. 

By putting "observers" at each apex of a cube of side 
i/2, where L is the simulation box size, we produce eight 
all-sky maps of resolution Ngide — 256 from simulation 
1. From simulation 2, we compute six 18° x 18° maps 
of resolution of Ngide = 2048, by putting an observer 
outside each face of the simulation box. Finally, from 
simulation 3, using the same setup as for simulation 2, we 
compute six 3° x3° maps oi Ngide = 8192. Figure [1] shows 
one map produced from each simulation; it should be 
noted that even the maps of patches of sky are computed 
using a spherical sky. 

Comparing our maps with those of pL5|, we note that 
the 18° map shown and their 7.2° map have similar fea- 
tures, but the former does not present as sharp linelike 
discontinuities as the latter. The resolution of our maps 
is effectively lower than that implied by the Healpix N^ide 
parameter used. This can be seen most directly from the 
power spectra (see Sec. |V|: Normally, one would expect 
(■max ^ iNside^ but the powcr in the maps falls around 
i « Nside, which shows they are over pixelized or, to put 
it another way, the effective resolution is about a third 
of that expected from the pixelization. For example, our 
18° maps have an effective resolution of about 5', com- 
pared to 1.7' expected from a map of Ngide = 2048. This 
explains the difference in features with the Fraisse et al. 
maps, which have a resolution of 0.42', an order of mag- 
nitude higher than ours. 



IV. MAPS 

We compute maps of CMB fiuctuations by following 
photon paths through the simulation boxes. The CMB 
temperature fluctuations are given by the following equa- 
tion, obtained by integrating the linear Boltzmann equa- 
tion for the Stokes parameter /: 

ST /"'" 

. (2) 

— ^e^'^hijUifij ] c?77 

where Ilf^ is the term that couples the Stokes parameter 
I to Q and U and is given in Fourier space by 

+{eiieij - e2ie2j)n'^+ + (eHe2j + e2ieij)Il'^^ 

(3) 

and all the other terms have their usual meaning. A full 
derivation of this formula is given in p3j . 



V. POWER SPECTRUM COMPUTATION 

For the all-sky maps, we decompose the temperature 
fiuctuations in the sky in spherical harmonics: ^ — 
"^2 (^imY™ and the angular power spectrum is estimated 
from 

= X! "^™"^m ■ (4) 

m 

For this purpose, we use the Healpix package [l^l- For 
the 3° p atches, we use the fiat-sky approximation (see 
e.g. [2^) which replaces the spherical harmonic trans- 
form with a 2-D Fourier transform: Ci ~ Cfc = |akp, 
where the modes are obtained from ^ = Qke*'^ ''. For 
the 18° patches, we have used both methods to compare 
their relative merit. To use a spherical transform on an 
incomplete sky, one must multiply the extracted spec- 
trum with a mode decoupling matrix [265: Cg — M^^lCi' , 
where 

Af«, = (2/ + 1) ^ iii^cr-^ (000)'' 
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FIG. 1: Temperature fluctuations produced by cosmic strings from simulation 1 (left), simulation 2 (middle), and simulation 3 
(right). 



However, the resulting spectra show strong oscillations 
due to the Gibbs effect that are not corrected by the 
procedure outlined above. Because of this, we have used 
a 2D Fourier transform even though the flat-sky approx- 
imation starts to break down at the lower end of the 
multipole range probed by these maps. In Fig. [21 we 
show the angular power spectrum from each of the three 
simulation sets, as well as their summation (i.e. concate- 
nating these consecutive but nonoverlapping time contri- 
butions). For the full sky, all multipoles are shown, while 
for the patches of sky, the spectra are binned. The errors 
on the individual parts of the spectrum are estimated by 
the variance between the maps. To add up the spectra in 
sections where two sets of maps contribute, we average 
the lower I part in bins of the same size as that of the 
higher I part and then add the two contributions. Also, 
to reduce the error, we binned the part of the spectrum 
in which only the all-sky maps contribute. We normalize 
the string spectra using Gji/c^ = 1 x 10~^. 

In addition. Fig. [3] illustrates the power spectrum from 
simulation 3 which begins before equal matter radia- 
tion and separates out the late-time contribution (dashed 
line), that is, as if the string simulation and Einstein- 
Boltzmann evolution were to start just before decou- 
pling ?7dcc- Note that maps by other groups have been 
generated by considering contributions only after decou- 
pling 77 > 77dec, e.g. ^5j. Here, we see starkly illus- 
trated the importance of the early contributions to the 
CMB anisotropy which arise from matter and radiation 
perturbations induced by strings before decoupling (pri- 
marily scalar modes). For / > 400, this strong early 
scalar component predominates over the late-time grav- 
itational string contribution (primarily vector modes). 
These results contrast markedly with the subdominant 
scalar power spectra obtained earlier using Nambu string 
simulations with full Einstein-Boltzmann evolution in a 
flat cold dark matter model [13] . However, this enhance- 
ment is due here to inclusion of the full matter-radiation 
transition with a significantly higher radiation era string 
density and a A concordance cosmology with a higher rel- 
ative baryon density at decoupling. In future work [28| . 



we shall compare this result with the spectrum obtained 
from the same simulations' UETCs. The overall shape 
for the spectrum computed from the 18° maps is also 
qualitatively consistent with the CMB power spectrum 
calculated using a simplified analytic formula [29i] . 



Comparing this result with that of [15| - the only other 
group to have recently extracted Ci from simulated maps 
- we see some significant differences attributable both to 
the inclusion of recombination physics in our work and 
to the different approximations used. Firstly, we note 
that Ref. fisj achieves much better statistics because it 
is easier to simulate purely gravitational effects in a flat- 
sky approximation at high resolution, generating a larger 
number of maps which results in smaller error bars. How- 
ever, for I < 500, we note that their power spectrum 
does not appear to turn over in the range I < 200-400, a 
causality cut-off due to the fact that these string modes 
are superhorizon at recombination (see e.g. Ref. 29]). 
In the range I 400-2000, as discussed above, there is 
an important contribution from early-time scalar modes 
at decoupling which is not included in Ref. [l^. It is 
reasonable to assume that late-time gravitational effects 
dominate the spectrum at I ^ 1000, and their result in 
this range is probably fairly accurate. However, in this 
region we note that the limited spatial resolution of our 
individual maps (resulting from propagation of photons 
through strings smoothed onto 256"^ grids) means that 
there is generically small-scale power missing at higher 
multipoles from the summed power spectrum in Fig. [2] 
For example, the gravitational contribution from strings 
in simulation 1 (full-sky) is absent beyond I « 100. For 
this reason and given the limited number of realizations, 
we have not yet endeavored to obtain a detailed con- 
straint using WMAP, though we expect it to be con- 
sistent with constraints in Ref. [s^]- By comparing the 
Sachs- Wolfe plateau to CMB observations, we find a nor- 
malization of the cosmic strings linear energy density to 
be Gii/(? = (1.45 ± 0.6) x 10^®, in agreement with our 
previous result [l3 |. 
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FIG. 2: On the left, we show the angular power spectrum of CMB temperature fluctuations produced by networks of cosmic 
strings The thin lines represent the spectra from simulations 1-3 (left to right) and the thick line is their sum. On the right, 
the total spectrum with its error bars is shown again as an "expectation area." 




FIG. 3: Power spectrum from simulation 3 using the full sim- 
ulation (solid line) and the late start evolution (dashed line) 
and their difference (dotted line). 



VI. PROBABILITY DISTRIBUTION 
FUNCTION 

One of the principal aspects of cosmic string induced 
CMB fluctuations is their intrinsic non-Gaussianities. In 
this section, we present preliminary results from our map 
realizations for the pixel temperature distribution, its 
skewness and kurtosis. In future work, we shall explore 
the efficiency of different techniques and tools, such as the 
bispectrum (sij . to better characterize the non-Gaussian 
signature from strings and thus infer their detectability. 

We compute the pixel probability distribution function 
of the rectangular map pixel grids computed to estimate 
the power spectrum. To compare our results, we gener- 
ate an ensemble of 1000 Gaussian maps with the same 
power spectrum as the string map. To do so, we use 
the fact that, in Fourier space, a Gaussian field will have 
phases that are uniformly distributed between and 2tt. 
Hence, for a given map, after fast Fourier transform, we 
randomize the phases and fast Fourier transform back to 
real space. In Fig. 21 we show the temperature distribu- 
tion for medium- and small-angle maps along with the 



Gaussian ensemble for comparison. These probability 
distribution functions are remarkably Gaussian, consis- 
tent with early string map results which demonstrated 
that, despite the distinct signature of individual strings, 
the central limit theorem prevailed after many strings 
contributed [11] ■ We note, however, that the higher res- 
olution 3° maps do appear to have a slightly increased 
level of non-Gaussianity. 

We also computed the skewness 
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and the kurtosis 
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(6) 



(7) 



of the string maps and compared it to that of all our 
Gaussian maps. Figure [5] shows the skewness and kurto- 
sis for all the maps as well as the Icr contour value for 
the ensemble average of Gaussian realizations. This in- 
dicates a marginal positive la skewness and kurtosis for 
the 3° maps. The absence of a strong skewness in the 
18° maps seems to be at variance with the significant 
negative skewness 71 = —0.24 found in the 7.2° maps 
from late-time gravitational effects in Ref. [l3| (see also 
the analytic estimates in [2^[3^[3l]). Two possible ex- 
planations are apparent. First, as mentioned above, the 
present 3D Einstein-Boltzmann simulations have an ef- 
fectively lower resolution than the 2D flat-sky approxima- 
tion maps, so they cannot probe as far into the wings of 
the distribution. We are missing the integrated effect of 
bispectrum triangles combining disparate large and small 
scales. Second, the analytic modeling in Ref. [2^ in- 
corporates a causality or correlation-length cut-off which 
prevents bispectrum and trispectrum contributions from 
superhorizon scales, effectively flattening their accumu- 
lated amplitude below I « 300 — 500. Without such a 
cut-off introduced through energy-momentum compen- 
sation, there may be a significant enhancement of the 
skewness. 



FIG. 4: Temperature distributions of the medium- (left) and small- (right) angle maps shown above. The solid lines indicate 
the la level of an ensemble of Gaussian maps 



The tentative indications of positive skewness and kur- 
tosis in the 3° maps, however, are consistent with phys- 
ical expectations. On scales above I > 300, we see that 
there is an important contribution from early scalars with 
strings generating wakelike objects in the matter distri- 
bution with a positive skewness and kurtosis [34|. This 
will lead to CMB anisotropics with similar properties. 
It appears that this positive skewness will compete with 
and confuse the negative skewness predicted for the late- 
time gravitational effects on intermediate scales. In con- 
trast, the kurtosis of these two contributions will sum 
positively, reinforcing the view that the trispectrum will 
prove a better discriminant of strings than the bispec- 
trum [29]. These preliminary results indicate that incor- 
porating recombination effects will be important in any 
effort to constrain non-Gaussian strings signatures with 
WMAP or Planck data. 



VII. CONCLUSION 

We have computed maps of cosmic string induced 
CMB fluctuations at various resolutions and we have ex- 



tracted their angular power spectra. We have demon- 
strated the importance of recombination effects for the 
power spectrum over a broad range of multipoles 200 < 
I < 2000. We have also shown that the resulting maps 
are remarkably Gaussian, though with potential devia- 
tions which are worthy of closer investigation as testable 
string signatures in the CMB. 
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